Precision spectroscopy of pionic Is states of Sn nuclei and evidence for partial 
restoration of chiral symmetry in the nuclear medium 
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Deeply bound Is states of 7r~ in ' ' Sn were preferentially observed using the Sn(d, 3 He) 
pion-transfer reaction under the recoil-free condition. The Is binding energies and widths were 
precisely determined, and were used to deduce the isovector parameter of the s-wave pion- nucleus 
potential to be b\ — —0.115 ± 0.007 m^ 1 . The observed enhancement of |6i| over the free ttN value 
(bi ee /bi = 0.78 ±0.05) indicates a reduction of the chiral order parameter, fn(p) 2 / fi ~ 0.64, at the 
normal nuclear density, p = po. 

PACS numbers: 36.10.Gv,13.75.Gx,14.40.Aq,25.45.Hi 
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Our work deals with characteristic properties of 
strongly interacting bound sytems, hadrons, consisting 
of light quarks (u and d). Their masses are nearly two 
orders of magnitude smaller (m u ~ 5MeV, md ~ 8MeV) 
than typical hadron masses of « 1 GeV. This extraordi- 
nary phenomenon is proposed to be produced by spon- 
taneous breaking of chiral symmetry for massless quarks 
subject to the strong interaction 0,0,0. It results in a 
ground state, the vacuum state, with a finite expectation 
value of quark- antiquark pairs, (qq)o ~ —(250 MeV) 3 Q. 
In such a scenario the hadrons are considered as quasi 
particle excitations of the condensate, (qq), separated by 
an energy gap of ~ 1 GeV. The lowest energy excitation 
modes of the condensate, so called Nambu-Goldstone 
bosons, are identified as pions. Their s-wavc interaction 
with nucleons is predicted to vanish in its isoscalar part, 
and determined by the pion decay constant, in its 
isovector part |5j, |fj which is consistent with recent ex- 
perimental values 0,13 ■ The f% is also the order param- 
eter of chiral symmetry breaking and directly connected 
to the magnitude of (qq) through the Gell-Mann-Oaks- 
Renner relation 

We examine this chiral symmetry scenario by im plant - 
ing a 7r~ in a nuclear medium of density p where 
a new vacuum state with a reduced condensate, (qq) p , is 
proposed to be created . This is the most basic case of 
a large effort to study the density and temperature de- 
pendence of the (qq) in relativistic heavy ion collisions. 

The density dependence of the quark condensate is ex- 
pressed in the leading order by the pion-nucleon cr-term 



(ctjv ~ 45 MeV) in the following form |Tl 



(gg)p 

(qq)a 
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which yields a reduction of about 0.65 for the nor- 
mal nuclear density, p = po = 0.17 fm~ 3 . Likewise, 
the pion decay constant in a medium (identified as the 

time component of the axial current) is reduced as [l^ . 

f * i \ 2 

p ~ 1 — otp, where the parameter a is predicted to 

be apo w 0.45 from the chiral dynamics This re- 

duction of the ratio is associated with the free and the 
in- medium isovector nN scattering length (b l { cc and &*, 
respectively) as 0] 
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Our program aims at measuring the isovector ttN 
interaction parameter in the pion-nucleus potential 
(b*(p)) [13 by studying deeply bound Is states of n~ 
in heavy N > Z nuclei B E E E 0. Such 
states were predicted to be pro duced as discrete states 
by nuclear reactions [2l|, H2, lH, |25| and to provide 
unique information on the s-wave interaction, whereas 
in most pionic atoms the repulsive s-wave interaction is 
nearly counterbalanced by the attractive p-wave inter- 
action (see, for instance, Batty et al. |2(j). In recent 
papers we showed that the density-dependent parameter 
b*(p) can be well represented by a constant parameter 
b\ |23, and developed a method to deduce the 61 pa- 
rameter El. The data of the Is n~ state in 205 Pb Hi, 
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combined with information on Is states of light pionic 



atoms, yielded 61 



-0.116 



+0.015 
-0.017 



A similar con- 



clusion was also obtained from global fits of pionic atom 
data by Friedman, bi = -0.122 ± 0.004 to" 1 [29| and 
— 0.113 ± 0.004 m™ 1 [3(J], although some ambiguity aris- 
ing from the p-wave part remains. In order to determine 
b\ more reliably and accurately, it is essential to perform 
high-precision spectroscopy on deeply bound Is ir~ states 
in heavy nuclei. 

In the present paper we report on systematic experi- 
mental studies of Is tt~ states in a series of Sn isotopes, 
which were produced with the Sn(d, 3 He) reactions. One 
of the advantages of using Sn isotopes is that we can 
produce the Is n~ states as the most dominant quasi- 
substitutional states, (ls) w - (3s)" 1 , because of the pres- 
ence of the 3s orbital near the Fermi surface, as theo- 
retically predicted [25j . Another merit is to make use 
of isotopes over a wide range of (N — Z)/A to test the 
isospin dependence |lfj| . 

We used a deuteron beam from the heavy-ion syn- 
chrotron SIS at GSI, Darmstadt, combined with the 
Fragment Separator as a high-resolution forward spec- 
trometer jlSj . The present experiment had many new 
features and improvements. We chose the exact recoil- 
less condition to suppress minor states other than the 
enhanced Is ir~ states with quasi-substitutional Is ir~ 
states, with a deuteron beam of a small momentum 
spread and an accurately measured energy of 503.388± 
0.100 MeV. Its intensity was 1.5 x 10 11 deuterons per 
spill. The beam was focused on a narrow spot with a 
width of about 1 mm, hitting a strip target with a width 
of 1.5 mm and a thickness of 20 mg/cm 2 . Enriched Sn 
isotopes with A = 116 (98.0 %), 120 (99.2 %) and 124 
(96.6 %) were used as targets. 

The strategy to arrive at high accuracy in determin- 
ing the binding energies and widths of the Is states 
was the following. First, we used the non-pionic 
A Sn(d, 3 He) A " 1 In(g.s.) (A = 112,116,120,124) reac- 
tions to calibrate the spectrometer momentum (with re- 
spect to the incident beam momentum) and dispersion 
to achieve an accuracy of 6 keV. The overall energy res- 
olution, which was crucial for measuring the widths of 
the Is 7T~ states, was also determined from these spec- 
tra to be 394 ± 33 keV (FWHM). Secondly, we put a 
thin mylar layer on the down-stream surface of each Sn 
target to measure the 3 He peak from the p(d, 3 He)7r° re- 
action as a built-in indicator of absolute calibration. We 
observed spectra, cPa/dE/dtl, on mylar-covered 116 Sn, 
120 Sn, 124 Sn targets as function of the 3 He kinetic en- 
ergy, as shown in Fig. ^ 

The skewed peaks at around 371 MeV arise from 
p(d, 3 He)7r°. Their kinematically broadened shapes were 
well explained by a simulation, which took into account 
the beam-energy profile at the thin mylar layer, the an- 
gular distribution, the spectrometer acceptance (A9 X rj 
±15 mr and A9 y ~ ±10 mr) and the kinematical energy 
shifts of the 3 He particles. It turned out that the energy 
at half maximum of the higher-energy edge of the peak 



(shown as vertical dotted lines in Fig. represents the 
kinetic energy of 3 He in the ideal case of 0-degree emis- 
sion from the p(d, 3 He)7r° reaction, and serves as an ab- 
solute calibration point with an accuracy of ± 7 keV. We 
adopted this "edge" calibration for our final results, since 
this long term calibration was pertinent to all runs of the 
three targets. It was found to agree within 13 keV with 
the short-term calibration using the Sn(d, 3 He)In(g.s.) re- 
actions. 
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FIG. 1: Double differential cross sections versus the 3 He 
kinetic energy of the 124 ' 120,116 Sn(d, 3 He) reactions measured 
at the incident deuteron energy of 503.388 MeV. The scales 
of the 7r~ binding energies in 123 . 119 - 115 g n are a ] so indicated. 

In each spectrum of Fig.^one distinct peak at around 
365 MeV was observed, which was assigned to a domi- 
nant configuration of (ls) w (3s)~ 1 . The overall spectrum 
shapes for the three Sn targets were found to be in good 
agreement with the predicted ones [25| . The spectra were 
decomposed according to the theoretical prescription of 
Ref. |25j |. from which we could precisely determine the 
Is binding energies (.Bis) and widths (As)- The results 
are summarized in Tableland in Fig. [21 The quoted er- 
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TABLE I: Observed binding energies (Bis) and widths (As) of the Is ty states in ' ' Sn isotopes. 

Igoto e B , McV1 ABis [MeV] , McV1 A As [McV] 

s stat. syst. total 8 stat. syst. total 

"""^Sn 31)06 ±0.021 ±0.012 ±0.024 0441 ±0.068 ±0.054 ±0.087 

119 Sn 3.820 ±0.013 ±0.012 ±0.018 0.326 ±0.047 ±0.065 ±0.080 

123 Sn 3.744 ±0.013 ±0.012 ±0.018 0.341 ±0.036 ±0.063 ±0.072 



rors are statistical, systematic (arising from the calibra- 
tion) and combined errors. In this study isotope shifts 
of deeply bound pionic Is states in heavy nuclei are seen 
for the first time. 
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FIG. 2: Experimental values of {Bis, As} (elliptic circles 
for statistical errors; horizontal and vertical bars for total 
errors) for tt~ in 115 . 119 > 123 Sn, shown together with theoretical 
relations with {fei, ImA)} for each isotope. 

The most serious problem in our analysis is the rela- 
tively poor knowledge concerning the neutron distribu- 
tion p n (r) in Sn isotopes, whereas the proton distribu- 
tion p p (r) is well known [3l|. We take the following form 
for the difference between the neutron and proton rms 
radii, which was derived based on experimental data of 
antiprotonic atoms of Sn isotopes |33 | as well as of many 
other nuclei Ar np = (1.01 ± 0.15)(iV -Z)/A+ (-0.04 ± 
0.03) fm. In converting Ar np into the diffuseness parame- 
ters (dp, o„) and the half-density radii (c p , c n ) in the two- 
parameter Fermi distribution we took half way between 
the "halo" type {c r > = c n ,a p < a n ) and the "skin" type 
(c p < c n ,a p = a n )|32], in accordance with a proton elas- 
tic scattering result in Pb isotopes (see Ref. :28|). The 
adopted parameters {c p , a p , c„, a n } in fm are: {5.405, 
0.523, 5.473, 0.552} for 115 Sn, {5.442, 0.523, 5.532, 0.560} 
for 119 Sn, and {5.484, 0.523, 5.594, 0.569} for 123 Sn. 

The influence of Ar np on the extraction of optical po- 
tential parameters in global fits of pionic atom data was 
recently discussed in detail by Friedman and Gal [3^ 



based on old theoretical values of Ar np [3_4J, which are 
considerably different from the empirical ones j^we use. 
These values are supported by a new proton-scattering 
result on Sn isotopes [35| and also by a new RHB calcu- 
lations |36j. In our view the use, in our analysis, of the 
experimental data available now resolves the concerns of 
Ref. [3j| with regard to our results. 

We made simultaneous fitting of B\ s and As of the 
three Sn isotopes together with those of symmetric light 
nuclei, leaving bo, bi, Rei?o and IuiBq as free parame- 
ters. Since the Is binding energy is insensitive to the 
p-wave potential, we could safely fix the p-wave param- 
eters to the known values from pionic atom data [26| . 
Hereafter, the units of m^ 1 for bo,bi and of m^ 4 for 
B will be omitted. We chose 16 0, 20 Ne and 28 Si for 
which the nuclear parameters are well known and the 
condition p p (r) — p n (r) is clearly fulfilled. Thus the ob- 
tained values are: b = -0.0233±0.0038, b x = -0.1149± 
0.0074, ReB = -0.019 ± 0.017, ImBo = 0.0472 ± 0.0013. 
The errors include both statistical and systematic uncer- 
tainties. In Fig. |2| we show the obtained relations of {£>i, 
Im£>o} mapped on {B\ s , As} for each Sn isotope. The 
intersections of the tilted axes are chosen to the best-fit 
values of {b\, IirBq}. We also show likelihood contours 
in the plane of {&i, lmB } in Fig. 

ap0 0.5 0.4 0.3 0.2 0.1 0.0 



bi" ee /bi 0.7 08 0.9 1.0 




FIG. 3: Likelihood contours in the {6i, ImA)} plane from 
the simultaneous fitting of {-Bis, As} of the Is pionic states 
in the three Sn isotopes and three light symmetric nuclei. The 
previous 205 Pb data reanalyzed with ImA) (= 0.046) is shown 
for comparison. 
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The smallness of ReB (-0.019 ± 0.017) supports the 
value (— 0.038±0.025) claimed from the requirement that 
the effective bo consists of the free irN value and the 
double-scattering term j2^. The combined isoscalar pa- 
rameter, b + 0.215 ReBo = -0.0274 ± 0.0002, is in good 
agreement with that obtained in the analysis of the 205 Pb 
[28| . The imaginary part, lmB = 0.0472 ± 0.0013, is 
consistent with the global-fit value of 0.055 ± 0.003 by 
Batty et al. [2(j and Friedman [2^, [3(j , considering that 
they included the angle-transformation (AT) term, which 
causes an appreciable decrease in the width [23 . In fact, 
the best-fit value in our analysis with the AT term in- 
cluded is: Imi? = 0.058 ± 0.003. These isoscalar pa- 
rameters are in fine agreement with those determined by 
using only the Is states of light symmetric pionic atoms. 
The addition of the Sn data consequently yields a precise 
value for b\. 

Figure EH also shows how the best-fit values move, if 
we choose the "halo" type (open square) or the "skin" 
type (open diamond) for p n {r) — p p (r). Furthermore, we 
find that the possible change of b u w hen we allowed an 
uncertainty of ±0.04 fm in Ar np |3j|, would be around 
±0.009, as indicated by two crosses in the figure. 

The magnitude of the observed \b\\ is significantly en- 
hanced over the free irN value, which translates into a 
reduction of/; 2 (Eq.EJ as 

R = b^ ee /h = 0.78 ± 0.05 (3) 



b f r/bi(p e )^f:(p e ) 2 /f^i-a Pe , (4) 



where we made use of the fact [53, that the solu- 
tion with a local-density-dependent parameter, b\{p) = 
bf ec /(l — ap(rj), is equivalent to that using a correspond- 
ing constant parameter b\ = 6^ rcc /(l — ap e ) with an ef- 
fective density p e ~ 0.6po- 

The above value hence implies that the chiral order 
parameter, f*(p) 2 , would be reduced by a factor of 
« 0.64 at the normal nuclear density p = p - Using 
Eq. (J2Jl in the analysis, we obtain an experimental value 
of apo ~ 0.36, as shown in Fig. 03 which is close to the 
value 0.45 predicted by chiral perturbation theory (l3j |. 
If a theoretical value, ps + 3 MeV (averaged over 
7r + and ir~ |l3jp. is inserted into an in-medium Gell- 
Mann-Oakes-Renner relation 0,0], (qq) Po /(qq)o will be 
(m^/m^) 2 x (1 — apo) ~ 0.67, which is in good agreement 
with the value of 0.65, as predicted by using Eq.JQl. We 
have thus found clear evidence for the partial restoration 
of chiral symmetry, probed by well-defined pionic states 
in a well- defined nuclear density. 
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